The Origin of Fermi Arcs in Cuprate Pseudogap States and Strong Constraints on Viable Theories 

of High-Temperature Superconductivity 
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A full Fermi surface exists in underdoped high-temperature superconductors if the temperature T lies above 
the pseudogap temperature T*. Below T* only arcs of Fermi surface survive, scaling with T /T* as T — > 0, with 
T* displaying strong doping dependence. There is no accepted explanation for this behavior. We show that 
generalizing the BCS theory of normal superconductivity to include rf-wave pairs and antiferromagnetism leads 
to the origin and doping dependence of the T* scale, and a quantitative description of Fermi arcs. These results 
place strong constraints on viable theories of high-temperature superconductivity. 



High-temperature superconductors were discovered more 
than 20 years ago [1], but there is no uniformly-accepted ex- 
planation for their properties |2j]. They exhibit pseudogaps 
(PG) at lower hole doping: partial energy gaps occurring be- 
low a temperature T* but above the superconducting critical 
temperature T c H , with T* and T c exhibiting strong and op- 
posite doping dependence for low hole-doping. It is widely 
believed that understanding the PG states and the T* scale is 
central to understanding the high-T c mechanism |2j,|4j]. 

For normal superconductors — described well by Bardeen- 
Cooper-Schrieffer (BCS) theory [5] — the Fermi surface 
(highest occupied fermion energy level) is key to understand- 
ing superconductivity. In optimally-doped and overdoped 
high-r c compounds a "normal" Fermi surface exists, and BCS 
theory utilizing af-wave singlet hole pairs seems applicable 
In underdoped cuprates this picture fails: PG states (ly- 
ing between T* and T c ) have anomalous Fermi surfaces, with 
spectral strength pushed away from the expected Fermi en- 
ergy by partial energy gaps. This is termed "gapping out", or 
(more loosely) "loss", of the Fermi surface. 

Angle-resolved photoemission spectroscopy (ARPES) [6, 
0] probes electronic properties Ja Si of states in high- 
temperature superconductors. ARPES data suggest a de- 
creased state density near the Fermi energy for temperature 
T < T* that is anisotropic in momentum, with strong T- 
dependence. A full Fermi surface is observed for T > T*; 
as the temperature decreases below T*, only arcs centered on 
the J-wave nodal lines (diagonals in momentum space) sur- 
vive |H, OH [H, scaling as T /T* to zero length for T -> 
lUlll . Thus gapping is anisotropic in momentum space, with 
an ungapped Fermi surface surviving only in the form of 
temperature-dependent Fermi arcs. 

Earlier theoretical work on Fermi arcs has employed a va- 
riety of approaches, including precursor pairs II12[ 1 1 311 - the 
perturbative renormalization group [ 14], high-temperature ex- 
pansions in the t-J model |l5f] , RVB models with strong gauge 
coupling 01611 . and time-reversal violating phases 11711 . but 
there is no agreed-upon explanation of Fermi arcs. The is- 
sue received renewed focus because of results by Kanigel, et 
al that place the strongest constraints yet on the nature 
of Fermi arcs. We report here a quantitative description of the 



Fermi-arc data in Ref . [11], with a theory that also consistently 
predicts the pseudogap temperature T* at the heart of the 
Fermi-arc scaling behavior. Thus, we offer a comprehensive 
solution to one of the biggest mysteries in high-temperature 
superconductivity: the origin of the pseudgap and of Fermi 
arcs in underdoped cuprates. More importantly, our results 
imply that only high-T c models with a strongly-constrained 
set of properties can be consistent with the detailed attributes 
of Fermi arcs and the scale T*. 

In earlier work we proposed an SU(4) model for the ground 
state properties of cuprate systems. Generalized SU(4) co- 
herent states are the simplest Hartree-Fock-Bogoliubov vari- 
ational solutions that incorporate antiferromagnetism (AF) 
competing with ti-wave sup erconductivity on a lattice with no 
double occupancy Il8l Il9l |20j Bill . They generalize BCS to 
incorporate AF self-consistently. If AF correlations vanish, 
the SU(4) gap equations reduce to the BCS gap equations; if 
instead singlet pairing vanishes, the SU(4) gap equations de- 
scribe an AF spin system with Neel order; for finite AF and 
SC correlations, the SU(4) gap equations describe the evolu- 
tion with doping (P) and temperature (T) of c/-wave Cooper 
pairs interacting with strong AF correlations. 

SU(4) states forbid double occupancy by symmetry, not by 
projection 111 SB . The exact zero-T ground state at half filling 
is an AF Mott insulator that evolves rapidly with hole doping 
into a superconductor with strong AF correlations in the un- 
derdoped region, and finally into a singlet, c/-wave supercon- 
ductor beyond a critical hole doping P q ~ 0.16-0.19 Q. At 
finite T for P > P q , pairing is reduced by thermal fluctuations 
and the pairing gap vanishes at T c . In the underdoped region 
(P < P q ) there are quantum fluctuations associated with AF 
correlations in addition to thermal fluctuations at finite T. The 
interplay between AF, pairing, and thermal fluctuations pro- 
duces PG states above T c in which the fermionic pairs interact 
through AF correlations. Thus the SU(4) model identifies the 
energy gap A q opened by the AF correlations as the pseudo- 
gap, while the (singlet and triplet fermion) pairs in these states 
may be viewed as of- wave preformed pairs. 

The SU(4) model can reproduce both the pseudogap and 
pairing gap with doping dependence quantitatively in agree- 
ment with data (see Ref. 112 110 . The SC transition tempera- 
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FIG. 1 : (Color online) SU(4) cuprate phase diagram compared with 
data. Strengths of the AF and singlet pairing correlations were de- 
termined by global fits to cuprate data I20L 12111 . The PG temperature 
is T* and the SC transition temperature is T c . The AF correlations 
vanish, leaving a pure singlet rf-wave condensate, above the critical 
doping P q . Dominant correlations in each region are indicted by italic 
labels. Data from Ref. 1 23] (arrows indicate lower limits). 



tures T c and PG transition temperatures T* are also well repro- 
duced, as illustrated in Fig.Q] This is our first significant find- 
ing: SU(4) coherent states give descriptions of pairing gaps 
and pseudogaps, and their transition temperatures T c and T*, 
that are in quantitative agreement with cuprate data. 
The 15 generators {^} of SU(4) are M = ~{n - £2) and 

£> f = L s gk c lf-kl 4 = L hk4+q/-k,j 

k k (1 ) 

Qi) = L C l+qfk,j Sij = £ c lj c k,j - jOfy 
k k 

and hermitian conjugates, where |, j, i, and j are spin indices, 
q = {%,%,%), the maximum number of doped holes or par- 
ticles that can form coherent pairs (assuming the half-filled 
normal state as vacuum) is £2, and s g k is the algebraic sign of 
the d-wave pair formfactor 



g(k) = g(k x ,k y ) = cosk x 



cosk y . 



(2) 



In these expressions M is the charge with n = Y, k n k tne elec- 
tron number operator, and njj, Qij, and Sij are tensor com- 
ponents of the triplet pair, staggered magnetization, and spin 
operators K, Q and S, respectively. 

These generators are summed over momentum k. They 
are appropriate for describing data that are not momentum- 
selected (for example, Fig. [TJ. However, ARPES Fermi-arc 
data exhibit explicit dependence on k. The SU(4) formalism 
may be extended to deal with this case by viewing the indi- 
vidual k components of the operators defined in Eq. ([]]) as the 
symmetry generators Sf (jfc): 

} = {D + (k) ,D(k), jft (k) , f (*), Q(k) , S(k) , n k } 



where, for example, the singlet pair generator D T in Eq. 
(Q3 is related to the it-dependent generators D^(k) by D' = 
L*>o^ (k), an d k > means either k x > or k y > 0. 

Instead of the global SU(4) symmetry generated by Eq. |[T), 
the symmetry now is a direct product of fc-dependent SU(4) 
groups, riifc>o ( X ) SU(4)t. The corresponding Hamiltonian is 



H 



E {Xkk'Q(k)-Q(k') 

k>0 k,k'>0 

+G^(k)D(e)+a^nHk) ■ %{k')}, 



(3) 



where £^ and are single-particle energies and occupation 
numbers, respectively, and the interaction strengths are 



Xkk'=x\(k)g{k')\ 



G%=G®8 k 8« 



('=0,1), 



where 8^ = 8(9)\go(k)\. The pair formfactor is expressed in 
terms of the magnitude of the hole momentum k and its az- 
imuthal angle 9: g(k) = g(k x ,k y ) = g(k,9), where |go(*)| is 
\g(k, 0)| at the antinodes (9 = 0,71/2), where it is maximal, 
and 8(9) = 1 except in a very narrow region near the nodal 
point, where it drops rapidly to zero at the node. 

We shall term this ^-dependent formalism the SU(4)j 
model. All results of the original ^-independent SU(4) model 
of Refs. d [H [U EH (including those of Fig. Q} are re- 
covered as a special case of the SU(4)# model for observables 
that are dominated by contributions from near the Fermi sur- 
face (k = kf) and averaged over k directions. However, general 
solutions of the SU(4)^ model give new ^-anisotropic proper- 
ties. Of direct relevance to the present discussion is that the 
temperature for the PG closure T*(k) becomes anisotropic in 
k; specifically, we derive in the SU(4)^ model 



T*(k) = T*(k h 9) = T* \g(k t ,9)/go(k f 



(4) 



where T* is the gap closure temperature measured by ARPES 
along the antinodal { 9 = 0, 7t/2) direction, which is the max- 
imum possible value of T* and should be somewhat larger 
than T* inferred from experiments that average over Jfc. See 
Ref. lr22Tl for details. Equation (0|i defines the full temperature 
and doping dependence of Fermi arcs. The ultimate physical 
reason for this result is that the singlet and triplet pairs inter- 
acting by AF interactions in the SU(4) PG state each carry a 
g(k) formfactor, which introduces a k dependence in the ef- 
fective AF coupling and thus in T*. 

In Fig. |2a) we show the behavior of g(k) as a function 
of (k x ,k y ). The components (k x ,k y ) or (kf,9) are constrained 
by the Fermi surface (P = k 2 ). Assuming an isotropic hole 
surface [dashed gray lines in Fig. [2 a)], we have 



(% - k x 



■ ky) 2 =k}=2n{\+P) 



(5) 



(we consider a more realistic Fermi surface below). For a 
given doping P (thus kf) and temperature T, with T*(k) = T 
by virtue of Eq. (01 under the constraint ©, we obtain the 
angles 9\ and 02 at which the PG closes (the angle 9 is defined 
in the inset to Fig. [2b)), and the length of the surviving Fermi 
arcis/U0 2 -0i| =k { A6. 
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FIG. 2: (Color online) The ^-anisotropic factor, pseudogap correlation energy, and graphical Fermi-arc solution, (a) The ^-anisotropic factor 
\g(k)\ (blue lines) and contours of equal hole doping P (dashed gray lines), in the first Brillouin zone, (b) Correlation energy ay(6) for 
unit a = kT* evaluated along the Fermi surface [curves of constant P in part (a)] as a function of the angle 6 (defined in inset). Curves 
corresponding to doping P = 0-0.3 lie almost on top of each other, indicating that y(0) is largely insensitive to doping, (c) Graphical solution 
of the Fermi-arc problem. The curve defines the PG correlation energy and horizontal lines correspond to constant thermal energy scales k B T 
associated with a given temperature T. Their intersections (black dots) represent points in the momentum space where the pseudogap is just 
closed by thermal fluctuations; these bracket arcs AO of surviving Fermi surface. 



This result may be interpreted graphically. In Fig. |2jb) we 
show y{9) = \g(kf,8)/go(kf)\ versus angle 9 along different 
Fermi surfaces [dashed gray lines in Fig. EJa)]. We see that 
7(0) is almost independent of doping P. Thus, solving Eq. 
with the Fermi surface constraint (O is equivalent to solving 



k B T*(k) = ay{9) 



a = k B T 



(6) 



with k B the Boltzmann constant. The PG correlation energy 
ay(9) depends on the k direction [Fig.|2c)]. The pseudo- 
gap closes when the thermal excitation energy k B T is compa- 
rable to the PG correlation energy. The intersections of hor- 
izontal lines of fixed k B T with the correlation energy curve 
[black dots in Fig.|2|c)] define Fermi-arc solutions 9\ and 62- 
Outside those points (solid blue curve), the correlation energy 
is larger than the energy of thermal fluctuations (shaded re- 
gion), the PG opens, and the Fermi surface is gapped. Inside 
these points (dotted red curve), the thermal energy exceeds 
the correlation energy, the PG closes, and the Fermi surface 
exists in an arc AO between the dots. When T > T*, the PG is 
closed in all directions and there is a full Fermi surface since 
k B T > a, and a is the maximum PG correlation energy. 

The arc solution exemplified graphically in Figs. |2jb)-(c), 
or algebraically in Eq. (O, is our second significant find- 
ing: For given temperature T, the anisotropy of y(9) parti- 
tions the £-space uniquely into regions having a Fermi surface 
[T > T*(k t , 9)] and those that do not [T < T*{k f , 9)]. As Fig. 
0c) suggests, arc lengths decrease with decreasing T at fixed 
doping, with a full Fermi surface at T = T*, but only the nodal 
points at T = 0. Doping dependence enters primarily through 
the maximum PG temperature T* and the Fermi momentum 
kf, the temperature dependence enters through T/T*. For 
fixed T, arcs shrink toward the nodal points with decreased 
doping because T* in T/T* increases at smaller doping (Fig. 



Q}. However, if Fermi-arc length is measured by the fraction 
A9 1 (k/2) and the temperature is scaled by T*, the weak dop- 
ing dependence of y(9) ensures that A9/(n/2) versus T/T* 
is almost independent of doping and compound. 

Kanigel, et al 111 111 report fractional arc lengths versus T/T* 
for Bi2212 that we plot in Fig. [3] Our theoretical solution for 
the fractional arc lengths [obtained from Fig. [2jc) or Eq. (0]i] 
is the solid line in Fig. [3] Agreement between data and the- 
ory is remarkable, given that the theoretical curve has no ad- 
justable parameters (it is determined completely by the param- 
eters fixed previously in Fig. |T|) and that we predict the scale 
T* implicit in the data with the same theory (Fig.[T). Note that 
the different behavior at high and low temperatures in Fig. [3] 
(rapid drop in arc length for decreasing T ~ T*, shifting to 
approximately linear decrease to zero arc-length at T = 0), is 
explained entirely by geometry in Fig.|2c). 

Realistic cuprate Fermi surfaces are flat near zone bound- 
aries. The dashed curve in Fig. [3] uses the flatter Fermi sur- 
face shown inset lower right. The similarity of dashed and 
solid curves indicates that arc solutions depend only weakly 
on differing curvature in nodal and antinodal regions. Abso- 
lute arc lengths depend on doping and temperature. However, 
the scaled arc lengths of Fig. [3] are near-universal functions 
only of the ratio T/T*, largely independent of compound, 
doping, and Fermi surface details, as data suggest. 

Figure[3]illustrates our third significant finding: agreement 
between theory and data with no parameter adjustment sug- 
gests that our SU(4) formalism outlines a complete solution 
of the Fermi-arc mystery. We now argue that this has impor- 
tant implications for understanding pseudogap behavior. 



The data in Fig.[3]have been interpreted [ 1 1] as representing 
a rapid drop from a full Fermi surface at T = T*, destroying 
the antinodal Fermi surface at essentially constant tempera- 
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FIG. 3: (Color online) Fermi arc length versus temperature. Exper- 
imental arc length is displayed as a percentage of full Fermi surface 
length vs. T/T* for underdoped Bi2212 [11]; the solid curve is our 
prediction for an isotropic Fermi surface (inset upper left); the dashed 
curve assumes a Fermi surface with flat antinodal segments (inset 
lower right). No parameters were adjusted to these data in either 
calculation and the curves are almost independent of doping. 



ture, followed by a linear decrease of Fermi-arc length with 
decreasing T/T* on the near-nodal region of the Fermi sur- 
face, extrapolating to zero arc length at the nodal points for 
T/T* — > 0. This suggests that nodal and antinodal Fermi sur- 
faces may be removed differently as T is lowered (see Ref. 
1 24]). For example, it has been argued [11] that abrupt re- 
moval of the antinodal surface at T ~ T* may be associated 
with a lattice vector connecting antinodal surfaces (Ref. lf24ll 
speculates that gapping of the antinodal surface may be as- 
sociated with charge ordering), while the linearly-varying re- 
moval of the nodal regions extrapolating to a nodal-point sur- 
face at T = may support a nodal liquid picture [25]. 

The present results invite simpler hypotheses. Our unified 
analysis suggests that Fig.[3]is consistent with removal of both 
nodal and antinodal surfaces by the same mechanism. There- 
fore, the qualitatively different variation of arc length with 
T/T* in nodal and antinodal regions is not sufficient to in- 
dicate that loss of Fermi surface proceeds by different mech- 
anisms in these regions. Nor do the data of Ref. Ill 111 im- 
ply unique support for nodal liquids. We have demonstrated 
specifically that the SU(4) model accounts quantitatively for 
Fermi arcs and 7**, without invoking nodal liquids. 

We have shown that a theory of Fermi arcs must make two 
correct predictions: (1) the scale T* and its doping depen- 
dence, and (2) that T*(k) oc y{Q). Any theory having a T*(k) 
consistent with Eq. (O can describe the scaled data of Fig. [3] 
if T* is taken from data. Hence, our fourth significant find- 
ing: scaled ARPES data (Fig. [3} test whether a pseudogap has 
a nodal structure similar to that of the superconducting state. 
But discriminating among different theories meeting this con- 



dition requires more: a quantitative, self-consistent descrip- 
tion of Fermi-surface loss and of the PG temperature scale T* 
and its doping dependence (Fig. [TJ. We conclude that only a 
highly-restricted set of models can be consistent with the ag- 
gregate properties of Fermi arcs. Finally, although our results 
indicate that an SU(4) mean field can account for Fermi arcs, 
some properties of the pseudogap state are expected to be in- 
fluenced significantly by SU(4) quantum fluctuations that we 
shall address in forthcoming papers. 

We thank Elbio Dagotto, Pengcheng Dai, and Takeshi 
Egami for extensive discussion, and Elbio Dagotto for call- 
ing our attention to an error in our original formulation. 
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